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Abstract 

This paper unravels the problem of an idealised pile-up of n infinite, equi-spaced walls of edge 
dislocations at equilibrium. We define a dimensionless parameter that depends on the geomet- 
ric, constitutive and loading parameters of the problem, and we identify five different scaling 
regimes corresponding to different values of that parameter for large n. For each of the cases we 
perform a rigorous micro-to-meso upscaling, and we obtain five expressions for the mesoscopic 
(continuum) internal stress. We recover some expressions for the internal stress that are already 
in use in the mechanical community, as well as some new models. The results in this paper 
offer a unifying approach to such models, since they can be viewed as the outcome of the same 
discrete dislocation setup, for different values of the dimensionless parameter (i.e., for different 
local dislocations arrangements). In addition, the rigorous nature of the upscaling removes the 
need for ad hoc assumptions. 
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1. Introduction 

Dislocations occupy a central position in discussions of the permanent deformation of metals 
because of their role as the main carriers of plastic deformation. Therefore it is necessary to 
incorporate their presence, or the main effect of their presence, in a plasticity model that aims 
for a predictive power. However, since the typical number of dislocations even in a small sample 
of the material is very high, formulating a model that keeps track of every single dislocation is 
out of reach except for very small-scale problems. This explains the interest in describing the 
collective behaviour of dislocations in terms of a continuum quantity: the dislocation density. 

The challenge in this scale transition consists in describing the time evolution of the disloca- 
tion density in a physically-driven way, by performing a rigorous upscaling from the dislocation 
scale to the dislocation density scale (also called meso-scale). This task has been mainly pursued 
phenomenologically or by means of a statistical mechanics approach and has produced a number 
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of competing models (e.g. [5, 6, 8, 11, 12, 13, 18, 21]). In the case of parallel, straight disloca- 
tions, the evolution of the dislocation density p(x, t) (here we assume, for simplicity, that all the 
dislocations have the same Burgers vector) is described in terms of a continuity equation of the 
following form: 

d,p + d x (pv) = 0, 

where the velocity v driving the evolution is of the form 

b 

v = - (o- - o- int ), (1) 

with b denoting the Burgers vector, B a linear drag coefficient, <x an externally applied shear 
stress (constant, for simplicity) and cr int = cr int (x,p, d x p) an internal stress accounting for the net 
effect of the interactions among dislocations. The different models available in the engineering 
literature differ in the expressions of the internal stress that they propose; the range of validity of 
the proposed theories is typically unclear, as well as the conditions under which one of them is to 
be preferred over another one. In fact most of the models have been derived phenomenologically, 
often starting from an ad hoc Ansatz, or assumptions have been made in the derivation that are 
not always justified. 

Here comes the essential difference of our approach: we obtain a continuum model for the 
dislocation density from a more fundamental discrete dislocation model using a rigorous math- 
ematical approach (r-convergence). As a consequence, we obtain an exact classification of lim- 
iting behaviour for the system that we study, which unifies existing, independently derived de- 
scriptions into a coherent framework, and identifies new regimes that have not yet been studied. 

In order to make each step in the derivation completely justified and rigorous, our starting 
point has to be an idealised dislocation configuration. More precisely, we consider the discrete 
model of an idealised pile-up of dislocations studied in [22] (see also [10, 15, 14, 2, 19]). This 
model describes the equilibrium positions of n dislocation walls under the influence of an applied 
stress cr that pushes the walls towards an impenetrable barrier; the barrier is modelled as a wall 
of pinned dislocations at x = (see Section 2 for the detailed description). 

The discrete equilibrium equations for the positions of the n walls can be written in the 
general form 

o J int -o- = 0, i = l,...,n, (2) 

where the discrete internal stress for the z'-th wall is the sum of the contributions due to the 
interactions with the other walls (see (3) for the detailed expression of the equations). Therefore 
the system (2) is nothing but v, = 0, where v, is the velocity of the z'-th wall. At this point one can 
intuitively imagine that passing to the (continuum) limit in the discrete equation (2) should give 
a continuum analog of (1) for v = 0, and consequently an expression for a continuum version of 
the internal stress. This is exactly the object we want to characterise. 

Our approach is different from the upscaling procedure followed in the quoted papers that 
considered the same pile-up configuration. In [15, 14] the convergence of stationary states is 
proved using formal methods, in a special case of our analysis. Moreover, while [2, 19] perform 
a two-step discrete-to-continuum upscaling by smearing out the dislocations first in the slip plane 
and then in the in-wall direction, we upscale in the two directions simultaneously. We stress that 
the continuum model we derive is in perfect agreement with our starting discrete model. We 
recall that in [2, 19], as the authors point out, the resulting continuum model has to be corrected 
to incorporate the missing interaction. 
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An interesting novelty in our result is that the mesoscopic internal stress we obtain does not 
depend on the dislocation density (and its gradient) only, but it also contains some more local 
information about the discrete arrangement of the dislocations, that the density alone would fail 
to capture. To be more precise, consider the two arrangements in Figure 1. They correspond to 
the same density, i.e., to the same number of dislocations per unit volume. But the dislocation 
patterns in the two cases are different and this will result in different expressions for the upscaled 
internal stress describing the overall interactions in the two cases. 




Figure 1 : Two different arrangements corresponding to the same density in the rectangle. 

This additional local information on the arrangement of the dislocations can be expressed 
in terms of an aspect ratio a, defined as the ratio between the typical distance Ax between two 
consecutive dislocations in the same row (slip plane) and the distance h between two consecutive 
dislocations in the same wall. In Figure 1 the first arrangement corresponds to a ~ 1, the second 
one to a <K 1 . 

The upscaling procedure we adopt for the micro-to-meso upscaling is based on a variational 
convergence called F-convergence, which is well-known in the mathematical community and has 
been successfully applied to a variety of problems in materials science, from fracture mechanics 
to homogenization, from magnetomechanics to dimension reduction. Our approach focusses on 
the discrete energy of the system of dislocations (minimisers of the energy are exactly the solu- 
tions to (2)). The discrete-to-continuum upscaling is done by letting the number of dislocations 
n become infinitely large. According to the different asymptotic behaviour of the aspect ratio 
a (i.e., according to the local distribution of the dislocations), five different expressions for the 
continuum energy can be derived (we refer to [9] for the details of the mathematical procedure); 
and, accordingly, five different expressions for the upscaled internal stress. 

The results we obtain show that the simplified discrete model taken as a starting point is 
not too simple. In fact the internal stresses resulting from our derivation are more general than 
many well-known models proposed in the engineering literature (see [7], [8] and [13]). More 
precisely, our result contains as a special case the internal stresses proposed in the quoted papers, 
explaining their range of validity. The comparison with previous models will be the subject of 
Section 9. 

The key advantage of our rigorous approach to upscaling compared to other methods is that it 
is exact. This means that once a discrete model is chosen (with its simplifications and limitations) 
the corresponding upscaled continuum model obtained by following our method is uniquely 
determined. 

Moreover, the fact that no simplifying assumptions are made during the derivation (but only 
at the level of the discrete model) explains why such a rich variety of continuum models can be 
obtained, even from a relatively simple model. This is not the case for other upscaling methods, 
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where the starting point can be more general, but a number of restrictions and simplifications are 
typically made in the derivation, affecting the generality of the upscaled models and their range 
of validity. 

The plan of the paper is as follows. In Section 2 we state the discrete problem and we com- 
ment on our special choice for the arrangement of the dislocations. An overview of the methods 
and of the results is given in Section 3, while the details of the upscaling are considered in Sec- 
tions 4-8. A comparison between our results and other well-known models in the engineering 
literature is the subject of Section 9, and a number of interesting comments on the the results of 
this paper are presented in Section 10. Section 1 1 contains a short conclusion. 

2. Statement of the problem 

The problem we are considering is the equilibrium of a system of n walls of straight edge 
dislocations (all with the same Burgers vector) under the action of an externally applied shear 
stress that pushes the walls towards an impenetrable barrier. Consecutive dislocations within the 
same wall are assumed to be equidistant, at distance h > 0, and we assume there are infinitely 
many dislocations in each wall. The impenetrable barrier is modeled as an infinite wall of pinned 
dislocations (the plane x = 0). We moreover assume that the walls are perfectly aligned (see 
Figure 2). Therefore the model is essentially one-dimensional, the only unknowns being the 
positions jci, . . . , x„ of the n dislocation walls. 1 




slip planes 



cr 



Figure 2: Pile-up of discrete dislocation walls. The walls are at positions i, > and the barrier is modelled as a wall of 
pinned dislocations at x = jfo = 0. 

While it is common in the literature to reduce to the case of straight dislocations and to 
consider a single slip system (rather than multiple slip or cross-slip), the assumption of having 
dislocations arranged in vertical walls may sound unnatural (even if used in many papers, see 



'From here on, tildes distinguish dimensional quantities from their non-dimensional counterparts; we will define 
non-dimensional positions ,r, below. 
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e.g. [2, 19, 15, 14, 22, 10]). It is clearly a simplification and we are aware of the fact that 
a complete analysis of a fully 3d-dislocation model would be preferable. On the other hand, 
this idealised configuration allows us to carry out a rigorous analysis, whereas more random 
configurations would require a statistical approach. Such a statistical approach, however, would 
typically require some additional assumptions in the derivation, and therefore it might well result 
in a limit model that is less general than the one rigorously derived from a more idealised setting. 

Moreover, although the discrete model we consider is highly idealised, it has a number of 
properties that make it both interesting and not unrealistic. The fact that multiple dislocations 
move along exactly the same slip plane is natural, because of the way they are generated from 
Frank-Read sources (e.g. [17, Sec. 8.6]). Also, although the assumption of an arrangement 
in equispaced vertical walls is clearly an idealisation, it is on the other hand not unrealistic 
since equispaced vertical walls are minimal-energy configurations. Walls of edge dislocations 
are locally stable, in the sense that if one of the dislocations deviates from its wall position, 
either horizontally or vertically, it experiences a restoring force from the other dislocations that 
pushes it back. Finally, the vertical organization in walls is also justified by correlation functions 
calculated from numerical simulations (e.g. [13]). 

For these reasons we believe that our simplified discrete system is not too simple and that it 
is worth exploring it. Moreover, the mesoscopic models we obtain are general enough to contain 
as a special case several well-known models in the engineering literature, including some models 
derived by statistical arguments. 

The equilibrium positions of the dislocation walls are obtained by solving the following equa- 
tions for every wall i = 1, . . . ,n: 

nGb Ixi- Xj 



V)M s r L )- tr= a (3) 



2h(l - v) 

where <p{s) = ^-^ — is the (globally repulsive) stress governing the wall-wall interactions, -cr is 
the constant applied shear stress, b is the length of the Burgers vector, v is the Poisson's ratio and 
G is the shear modulus. Setting 

= nGb 
2(1 -v)' 

we can rewrite (3) as 

jr " I x- — x \ 

*Z^hr i - or=0, i = h -' n - (4) 

j*i 

Solving the equations (4) numerically gives a vector x e W; the i-th component of the vector 
represents the equilibrium position of the i-th dislocation wall. Defining the discrete dislocation 
density as 

P d (xd := . \ , i=l,...,n, (5) 

X[ X[— i 

we can plot the discrete density as a function of the equilibrium positions, obtaining the plots in 
Figure 3 for different values of the applied stress cr. 

From Figure 3 we can clearly see that the solution highly depends on the values of the prob- 
lem parameters. In particular, for large values of the applied stress cr the peak of the dislocation 
density close to zero seems to be dominating, while for smaller values of the stress a rather linear 

5 




Figure 3: Plots of the discrete dislocation density defined as in (5) corresponding to solutions of (4) for different values 
of the applied stress cr. In all the plots K = 1, n = 150, and h = 10; in (a) cr = 40; in (b) a = 0.01; in (c) a = 0.0005. 

behaviour of the density is prominent at first, and then a decay. This suggests that depending on 
the problem parameters, some behaviour of the density seems to be dominant in the bulk (suffi- 
ciently far from the boundary of the pile-up region, where the occurrence of boundary layers is 
expected). 

This key observation is already present in the recent paper [10], where a careful numerical 
analysis of the behaviour of the dislocation density in a pile-up is performed, and three different 
parameter regimes are identified. In particular, the linear dependence of the dislocation density 
with respect to the positions of the dislocations was observed there. 

However, a discrete numerical approach has some drawbacks. First of all, it is computation- 
ally expensive when the number n of dislocation walls is large. But the main drawback is that it is 
too detailed: knowing the position of each dislocation is definitely not necessary, i.e., an average 
mesoscopic information in terms of e.g. a continuum dislocation density should be preferred. 

3. Discrete energy formulation, scalings, and main results 

Despite its simple formulation, computing the continuum limit of the equilibrium equa- 
tions (4) turns out to be a rather involved problem. The complexity of the problem is illustrated 
in Figure 3, where the optimal discrete dislocation density (obtained numerically) is plotted for 
different values of the problem parameters. 

The aim of this paper is to derive a continuum description from the discrete system by rigor- 
ous mathematical upscaling. We do this by exploiting the variational structure of the equilibrium 
equations (4); hence we write down the discrete energy whose minimisers satisfy the equilibrium 
equations (4). The upscaling procedure will then be carried out on the discrete energy rather than 
on the discrete equations, and will produce a continuum energy functional. The method used 
to perform the discrete-to-continuum upscaling is T-convergence and the details of the rigor- 
ous mathematical derivation can be found in [9]. The equilibrium equation associated with the 
upscaled limit functional is exactly the upscaled version of the discrete equations and hence char- 
acterises the mesoscopic internal stress as desired. The steps of this procedure are illustrated in 
Figure 4. 

We now give an outline of the method and the results. We will discuss these results and put 
them in perspective, in Sections 4-10. 
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Figure 4: Main steps of our upscaling procedure. 



3.1. The discrete energy 

We define the energy associated to the equations (4). For ieR" let E„(x) be defined as: 



„ n n 



!=0 j=0 



(6) 



i=0 



where V is the primitive of -<p which decays to zero at infinity, defined as 

V(s) := -scothns — \ log(2sinh7Ti), 

7T 7T Z 

and plotted in Figure 5. The equations (4) are equivalent to d^E n = 0. 



(7) 



The first term in the energy (6) penalises configurations where the dislocations are close to 
one another, since the density V blows up logarithmically in zero, while it favours configurations 
where dislocations are at large distance from one another, since V decays exponentially at infinity. 
The second term of the energy, on the contrary, penalises configurations where the dislocations 
are far from the obstacle at x = 0. Using the fact that V(s) = V(-s) we can rewrite the sum in (6) 
as a one-sided sum, as 



n n-k 



E n (x) = Kj]J]v 

k=l 7=0 



Xj+k Xj 



(8) 



The form (8) of the energy will be more convenient for our analysis. 



3.2. Reseating 

We convert (8) into a dimensionless form involving the smallest possible number of indepen- 
dent parameters. This procedure will provide the relevant energy scalings that we consider in the 
discrete-to-continuum derivation. 

We first non-dimensionalize the dislocation positions x, by a length scale £„, by introducing 
the new dimensionless variables x,, 



0.8 



0.7 




Figure 5: The interaction energy V. 



At this stage the length scale £„ is undetermined; we will use our freedom of choice of £„ below 
to capture the behaviour of the pile-up. This is akin to choosing the appropriate degree of mag- 
nification in an experimental setup. The choices that we make below will make l n of the same 
order as the pile-up length. 

Given a choice of {„, the aspect ratio a n is the ratio between the average dislocation distance 
Ax ~ — and the vertical spacing h between consecutive dislocations in the same wall: 

r 

a„ := — . 
nh 

We now rewrite the energy (8) in a dimensionless form, in terms of the new variable x: 

E " (X) ^ = ^hVn 2 V K^+* " *J>) + l Z *J- 

k=l j=0 j=0 



(9) 



Since E„ will be of order one, the dimensional energy will scale as n 2 crha n . 

The aspect ratio a„ is a local description of the average arrangement of the dislocations within 
the pile-up region, whereas ( n gives a more global information on the dislocation distribution. 
Therefore there are two ratios of length scales: £ n /h, which compares the total length of the pile- 
up with the vertical spacing h, and a„, which compares the spacing between the walls and the 
spacing within the walls. 

The possible relations between these two (clearly non-independent) parameters is illustrated 
in Figure 6. The region marked as inaccessible corresponds to the regime a„ > %, which is not 
compatible with the definition of a„ since it would require that the length of the pile-up region is 
smaller than the spacing between the walls. 

The scalings of a„ illustrated in Figure 6 and that we will describe are listed below: 



(1) Subcritical regime: a„ <K 

(2) First critical regime: a„ ~ 



Figure 6: Scalings (l)-(5) in the parameter space (^-. Q 'n) an <i some of the corresponding configurations. 



(3) Intermediate regime: 1 •« a„ <K 1; 

(4) Second critical regime: a„ ~ 1; t ^ » 1 

(5) Supercritical regime: »„ « 1. 

3.3. Determining the scaling by a trial configuration 

In the companion paper [9] we prove that, given certain choices of €„ for each of the cases 
(l)-(5) above, the rescaled energy E„ converges as n — > oo. Here we do not give the details of 
this proof; instead we only specify the choices of {„, but motivate, heuristically, why these are 
the right ones, and how the resulting limit expressions for the energy and the internal stress arise. 

A convenient method to determine the scaling of t n is by requiring that both terms of the 
energy (9) are of the same order. As it turns out, the correct scaling can be determined by only 
considering the special, uniformly spaced configuration x, := i/n. In this case, the second term 
in (9) is always of order one, since 




as n — ^ oo, 
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and therefore we will require that the choice of £ n makes the first term in (9) also of order one. 

By this analysis we are able to relate („ and a n to the parameters h, cr, K and the number n of 
the dislocation walls. In particular, we relate („ and a„ to the dimensionless parameter/^ defined 
as 



introduced in [10]. In mechanical terms, j3 n is strongly related to the forcing term cr and is small 
for large forcing and large for small forcing. 

Note that this implies that K, cr, and h may depend on n, and in many cases this makes sense. 
For instance, if the number of dislocations n increases because the density of Frank-Read sources 
increases, then one would expect that the average slip plane spacing h should simultaneously 
decrease with n. 

Performing this matching on the trial configuration of uniformly spaced dislocations, as we 
do in Sections 4—8 below, we find that both terms in (9) have the same scaling provided we make 
the following choices: 

(a) If /?„ <K 1 /n, then choose 



Since a n = /3 n , these three cases correspond to the First critical regime, the Intermediate 
regime, and the Second critical regime, respectively, 
(c) If /3„ » 1, then choose 



Then a n » 1, and this is the Supercritical regime. 

In this way we are able to determine in which regime the system is by considering only /?„, which 
is an explicit function of the constants, and independent of the solution. 

3.4. Internal stress 

The most convenient concept to connect the discrete (finite n) situation with its infinite-n 
limit is the density p d introduced in (5). As n becomes large, this discrete density approximates 
a continuous density p. For four of the five regimes we derive expressions for the internal stress 
in terms of this limiting, continuous density p: 




(10) 



(b) 
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(1) Subcritical regime: cr\ x l{x) - —r I log \x - y\ dyp(y) dy; (11) 

n 2 Jo 

(2) First critical regime (assuming f3„ = 1/n): <r ( ^ t (x) - -K I v( — - — ) dyp{y) dy; (12) 

J h 



Kh 



(3) Intermediate regime: cr' ;(x) = 3i/3; (13) 

(4) K 1 , / 1 \ 

(4) Second critical regime (assuming /?„ = 1): cr| nt '(x) = —— tp eS I — I c^p. (14) 
Here 

CO CO 

jfc=l k=\ 

and has the asymptotic behaviour (see Section 7) 

1 



as f -> 
3^-f 3 



-87rte 2m as t -> oo. 

Note that the argument of <p' eS , 1 /hp, is an aspect ratio, since it is the ratio of the spacing (x ; — jc,-_i) 
to /z. 

The fifth, Supercritical, regime, has no meaningful internal stress, since it is too degenerate 
(see Section 8). 

3.5. Pile-up behaviour 

For each of the four cases above, the equilibrium density p of a set of dislocation walls that 
is pushed against the obstacle at x = solves the equation 

cr int (x) - o- = 0. (15) 

1. For the Subcritical regime, equation (15) has been solved by Head and Louat [16]: 



, cr \l - x . , „ 2nK 

^ = ?V— • with*=— . 



2. For the First critical regime, equation (15) has no explicit solution, to our knowledge. 

3. In the Intermediate regime, the equilibrium density is given by 



pi x)= ¥ (t-x), with* = y — . 

4. In the Second critical regime, again we believe that equation (15) has no known explicit 
solution. 

5. Finally, in the Supercritical regime, the equilibrium density p is constant: 

, [T 1 if x< €, . , „ 1 , / 2K x 

p(x) = { with£ = — n z ho-ln( . 

10 otherwise, 2n ^nhcr' 

11 



In the three explicit cases above, the parameter t is exactly the pile-up length. Note that I scales 
the same as the scaling parameter {„, but the two differ by a numerical constant. In Figures 8-14 
the continuum pile-up densities are compared to their discrete counterparts. 

We now consider each of the five regimes separately in Sections 4-8; in every case we vali- 
date the continuum upscaled model by comparing its predictions with the solution of the discrete 
model. A comparison between the internal stresses cr* above (for k = 1 - 5) with other well- 
known models in the engineering literature will be the subject of Section 9. Section 10 contains 
a number of interesting comments on the the results of this paper. 



4. Subcritical regime a n « \ : derivation of the upscaled internal stress 

In this section we consider the extreme case a„ <K L , or equivalently, f <K 1, i.e., the case 
where the length l n of the pile-up region is much smaller than the in-wall dislocation spacing. In 
this case the in-plane interaction is much stronger than the in-wall interaction, i.e., the equivalent 
continuum formulation will not sense the walls and the result will correspond to the case of a 
single slip plane rather than infinite walls of dislocations. 

This scaling regime has been analysed numerically (in its discrete formulation) in the pa- 
per [10], and corresponds to the so-called close region in their terminology. 

4.1. Heuristics for the scaling of the discrete energy 

As described above, we will identify ( n and a„ and the corresponding rescaling of the energy 
in (9) by requiring that the two terms of the discrete energy E„ are of the same order and bounded 
when calculated for the uniformly-spaced configuration. 

Since the second term in (9) is of order one, the only condition we are left to impose is that 
the first term of the energy is of order one, i.e., 



K 



J]v(a n k)~L (16) 



ncrha n h 

Since for this rescaling na n «: 1 , the argument of V in ( 16) is small, as a n k < a n n «c 1 . Therefore 
we substitute for V in (16) its asymptotic expansion close to zero, 

1 - log 2ns 
V(s) f , fors>0. 

n z 

Using this approximation in (16) we find that we should require that 

1 - log(2na n k) 



— 7 

ncrha n f-j 



~ 1. (17) 



As it stands, this condition would give rise to a wrong choice of ( n and a n . This is because the 
expression between brackets is dominated by a constant term, and this constant term is irrelevant 
for the equilibrium equations (since it vanishes upon differentiation). We therefore introduce the 
renormalized energy density V n (t) := V(t) H — 1+ ° % ^ nna " ) in which we have subtracted this large 
constant. In terms of V n the bound (17) on the energy becomes 

1 K A /*> 



n 1 ncrha n . 
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Simple computations show that 



K 

ntr ha, 



n 0g \nj o-ha„ [n \n) 



K 

(rhan 



(18) 



since the term in square brackets is the Riemann sum for the integral J Q log(f) dt = 1 . Hence 
the bound on the energy reduces to the requirement that ~ 1. This bound provides the 
following expressions for the aspect ratio a n and for the length l n of the pile-up region in terms 
of the parameters cr, h, K and n: 



K 



Kn 



crh' " cr 
or, in terms of the dimensionless parameter /3„ defined in (10), 



(19) 



(20) 



We note that, by (20), the scaling regime a„ <K - can be equivalently formulated in terms of f3 n 
and corresponds to /3„ <sc From (9), using (19) and (20), we obtain the scaling of the discrete 
energy in the subcritical regime, which reads 



n n-k 



£<»(*) : = 1 J] % (n 2 ti(x j+k - xj)) + - j] *J- (21) 

*=1 7=0 " 7=0 

Note that V„ can be expressed in terms of /3 n as V n (t) = V(t) + ±^S^lMl_ 

4.2. Continuum limit: derivation of the internal stress 

In [9] we rigorously derived a continuum energy in terms of a continuum dislocation density 
p, starting from the discrete energy (21), in the limit n — > oo. Here we give a hint of the main 
idea of the proof, and we refer to [9] for the detailed derivation. 

Starting from the discrete energy (21) we define p„ := - £"=i 6 Xi , where 8 Xi is the Dirac delta 
function localised at x, (and is zero everywhere except at x,). The measure p„ describes the 
distribution of the walls and is approximately a dimensionless version of the discrete density p d 
introduced in (5). In terms of p„ we can rewrite the sums in (21) as integrals, namely 

E ( n\x)=-\ V n (n 2 l3 2 „(x-y))p n (x)p n (y)dxdy+ xp„(x)dx. (22) 

z Jo Jo Jo 

In the previous formula we wrote the integral on (0, oo) since the wall positions x, are the un- 
knowns we are solving for, which are in (0, oo). We note that in terms of p n the discrete energy 
looks like a continuum functional. By the definition of V„ and since n/3„ «c 1 in this regime, we 
have that V„ {n 2 p 2 sj ^ log Substituting this expression in (22) we have 

E n\ x )--^~i J J \og\x -y\p„(x)p„(y)dxdy + J xp„(x)dx. (23) 
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It turns out that for a large number of walls, i.e., as n — > oo, p n converges to a continuum density p 
and the discrete energy converges to the continuum functional defined as: 

£ (1) (P):=-^ J o J \og\x -y\p(x)p(y)dxdy + xp(x)dx. (24) 

This is the top arrow in Figure 4. 

According to the scheme in Figure 4, it now remains to compute the Euler-Lagrange equation 
associated with the continuum functional £ ,(1) . In this case, the continuum equilibrium equation 
is an integral equation of the following form: 



i r 

~^ Jo 



P(j) dy+l=Q, (25) 



\x-y\ 

for every x e (0, oo). This equation is to be interpreted as a Cauchy principal-value integral (see 
e.g. [20, Ch. 2]). Alternatively, we can rewrite the equation (25) in the compact form 

-^log*^p + l =0. (26) 



We recall that our starting point was a discrete system of the form <x mt - cr — 0. Therefore the 
continuum equation (26) can be interpreted as the dimensionless form of - cr = 0. Hence, 
the expression for the dimensionless internal stress we obtained is 



o™ = ~log*d x p. (27) 



After returning to dimensional variables we find (1 1). 



4.3. The Eshelby-Frank-Nabarro model (EFN): the single slip plane case 

Since in this regime the length of the pile-up region {„ is much smaller than the vertical 
spacing h, the in-wall interaction is much weaker than the in-plane interaction. Therefore it is 
natural to compare our model with the EFN model, formulated in the case of a single slip plane. 

In the EFN model the equilibrium positions of n dislocations (as opposed to walls) in a single 
slip plane (under an applied stress cr that pushes them towards the barrier) are described by 
equations similar to (3). In that case, the interaction potential is iff(t) := 1 /n 2 t. 

In Figure 7 we show the comparison between the EFN model and our model. We plot the 
discrete density p d corresponding to the minimisers x e W of the energy (21), for /3 n <K -, and 
of the EFN energy, where p d (x,) := l/(x, - x,_i). The discrete densities corresponding to the two 
models show perfect agreement. 

4.4. Comparison Discrete vs Continuum 

In this section we numerically compare the discrete density obtained by minimising the dis- 
crete energy (21) for /?„ <sc - and for large n with the solution of the continuum equation (26). 
The agreement is shown in Figure 8. 

We note that the continuum equation (26) can be solved by means of the Hilbert transform, 
and it has a closed-form solution, found in [16] (see also [24]). 
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Figure 7: Comparison between the optimal discrete density in the EFN model and in (21), in the subcritical scaling 
regime (1). Here n = 150 and /J„ = ^4= (<K lj. The stars and crosses completely coincide. 
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Figure 8: Comparison of the discrete and continuous pile-ups, for n = 150 and /?„ = 6/(n yfn). The continuum dimen- 
sionless density p minimizes the discrete dimensionless density p' 1 is calculated by minimizing e\, ' over all x,, and 
defining p d (Xj) := l/C*, - Xj-i). Both are normalized such that the total density equals 1. 
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5. First critical regime a n ~ 1 : derivation of the upscaled internal stress 



This regime corresponds to £ n ~ h, i.e., to the case where the length of the pile-up region is of 
the same order as the vertical dislocation spacing h and consequently, for large n, the horizontal 
wall-wall spacing is much smaller than h. Unlike the previous case, we cannot expect that the 
vertical interactions can be neglected and that the stress exerted by a wall is equivalent to the 
stress generated by a single dislocation — at least, not quantitatively. Qualitatively, though, the 
optimal dislocation density exhibits a sharp increase close to the obstacle, as in the previous case, 
and a fast decay at infinity. 

5.1. Heuristics for the scaling of the discrete energy 

To find the scaling of the discrete energy E n in (9) that guarantees that both terms are bounded 
and of the same order (and the corresponding expressions for a„ and £„) we use as a test con- 
figuration, as in the previous scaling, the uniform distribution x, ■ - l -. Since with this choice the 
linear term in the energy is of order one, the bound on the energy reduces to a bound on the first 
term only, i.e., proceeding as in (18), 



g n n-k g n 

i ~ - 1 - r - y y v (<*„*) — — 2 a n y v (<*„*) 

n L crha n *r-£ ncrhai. j—f 

k=l j=0 " L k=l 



(28) 



Since for this scaling na„ ~ 1 as n — > oo, the term in square brackets is of the same order as the 
Riemann sum for the integral J Q V(t)dt < oo. Therefore, the condition (28) reduces to 



1, (29) 



ncrha 2 , 

which leads to the following expressions for a„ and t n (by using also (10)): 



K "'' n/3 n h. (30) 

cr 



Note that by (30) we can directly reformulate the scaling regime as /3 n ~ K More precisely, this 
corresponds to /3 n = ^, for some c n ~ 1 (c n — > c as n — > oo). The rescaling of the energy (9) 
obtained in this case is therefore 



n n-k 



£ « 2) w : = % Z Z v ( Cn(x ^ - x $ + \ Z x i- (31) 

" t=i j=o n j=0 

5.2. Continuum limit: derivation of the internal stress 

We rewrite the discrete energy (31) in terms of the empirical measure p n - \ Zj=i an d 
we obtain 

E „ 2 \ x ) = TT | I V(c n (x - y))p n (x)p n (y)dxdy + I xp n (x)dx. (32) 
1 Jo Jo Jo 

For a large number of walls, i.e., as n — > oo, p n converges to a continuum density p, c„ — > c, and 
the energy E n 2 \x) converges to the continuum functional E (2> defined as: 

V(c(x-y))p(x)p(y)dxdy+ xp(x)dx. (33) 
o Jo Jo 
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The equilibrium dislocation density p is the solution of the Euler-Lagrange equation associ- 
ated with the functional £ (2) , i.e., the solution of the integral equation 

V(c(x-y))d y p(y)dy + l= (34) 

o 

for every x e (0, oo), or, equivalently, 

V c * 8 xP +1=0, (35) 

where V c {s) := cV{cs). This equation is again the mesoscopic equilibrium equation cnu-cr = in 
its non-dimensional form. Hence, the dimensionless internal stress obtained from this rescaling 
is 

o%=-V c *d x p. 

The dimensional version, for c — 1, is (12). 

We note that, as for the previous scaling regime, also in this case the continuum equilibrium 
equation is a singular integral equation. The main difference is that whereas in the rescaling 
case (1) we could approximate V with its limit behaviour near zero (namely its logarithmic be- 
haviour), here the complete energy density V enters the limit functional (33) and the equilibrium 
equation (35); the scaling constant c = lim„^ M n/3„ enters the expression as well. 



5.3. Comparison Discrete vs Continuum 

In this section we show the agreement between the solution of the upscaled continuum equa- 
tion (34) and the minimiser of the discrete energy (31), for a large number n of dislocation walls. 
The plot is shown in Figure 9. 
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Figure 9: Comparison of the discrete and continuous pile-ups for and £* 2 ', for n = 150/3,, = 5/n. 
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6. Intermediate regime 1 «. a n «. 1: derivation of the upscaled internal stress 



This scaling regime corresponds to the intermediate situation in which the pile-up length ( n is 
much larger than the vertical spacing h while the average in-plane distance between consecutive 
dislocations is, on the contrary, much smaller than h. This case corresponds to the remote region 
considered in [10]. 

6.1. Heuristics for the scaling of the discrete energy 

Using again the equispaced wall distribution x, = L as a test for the discrete energy E„ in (9) 
and proceeding as in the previous cases leads to the following requirement 



n 2 crha n 4-i " ncrhal 

k=l j=Q " 



" 2 V (""^ 



k=\ 



(36) 



since, by assumption, a n — » and na n — > oo, the term in the square brackets is a Riemann sum 
for the integral J Q V(t)dt = l/6n. Hence the condition (36) on the energy is equivalent to the 
following relation among the parameters: 

1, (37) 



ncrhal 



which is identical to (29) and therefore leads to the same expressions for a n 3) and 6^ as for a4 2) 
and 4, 2) in (30). In particular, also in this case a n 3) ~ f3„; in terms of /3„ the scaling regime is 
- <K /3„ «: 1 . This relation between 4? and the parameters h, cr, K and n found above was also 
found in [10]. 

Using (37), the rescaling of the energy (9) becomes 

~ n n-k . n 

E n 3) (x) := ^ J] Yj V { n Pn(xj+k - xi)) + - Y xj, (38) 

k=\ 7=0 7=0 

for \ «/3„ « 1. 

6.2. Continuum limit: derivation of the internal stress 

In terms of the distribution p n -\ X"=i S Xi , the energy (38) can be rewritten as 

E^\x)=-np n \ V(n/3 n (x-y)) Pn (x)p n (y)dxdy + x Pn (x)dx. (39) 

1 Jo Jo Jo 

We note that, since n/3„ — > oo, n/3„V(nfi n s) — > (JT^ v) <5o = (3;r) _1 5o; hence for n — > oo the energy 
£^ 3) converges to the continuum energy £ (3) defined as: 

£ (3) (p) := r" I P 2 (*V* + | *p(*)dx. (40) 
6?r Jo Jo 

The dislocation density p minimising £ (3) is the solution of the Euler-Lagrange equation associ- 
ated to the dimensionless functional £ (3) , i.e., 

^d x p +1=0. (41) 
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We notice that unlike the previous scaling regimes the continuum equilibrium equation in this 
case is local. Moreover, from (41) it is clear that the optimal dislocation density is linear. This 
was observed numerically in [10] and proved in [14] using formal methods. 

(3) 

Equation (41) is the mesoscopic equilibrium equation cr>' - cr — in its non-dimensional 
form. Hence, the dimensionless internal stress obtained from this rescaling is 

(3) * a 
a: i = o x p. 

6.3. Comparison Discrete vs Continuum 

In Figure 10 we compare numerically the discrete density corresponding to a minimiser of 
the rescaled energy E„ in (38) and the solution of the continuum equation (41). We notice that 
the agreement between the two densities is perfect in the bulk, sufficiently far from the boundary 
of the pile-up region, where the occurrence of boundary layers is expected. 




15 - 



10 




0.05 0.1 0.15 0.2 0.25 0.3 
Dislocation wall positions 



0.35 



0.4 0.45 



Figure 10: Comparison of the discrete and continuous pile-ups for E„ and for n = 150 and fj„ = 1/ -*Jn = 1/ V150 



6.4. Comparison with the Evers-Brekelmans-Geers continuum model. 

The internal stress cr. coincides, at least qualitatively, with the one proposed in [8] and 
derived phenomenologically from interactions among dislocations. In the quoted paper, how- 
ever, the expression of the internal stress contains a length scale R (not explicitly determined), 
representing the spatial reach of the dislocations interactions. By matching the stress cn^ in its 
dimensional form (13) with the one derived in [8] we can determine the length scale R, namely 
R ~ h. Therefore R depends on the in- wall spacing h only. 
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7. Second critical regime a n ~ 1: derivation of the upscaled internal stress 



This scaling regime corresponds to configurations where the in-plane and the in-wall dis- 
tances between consecutive dislocations are of the same order. As a result, the length of the 
pile-up region is much larger than h. 

The scaling a„ ~ 1 is a special situation as compared to the cases (l)-(3) treated so far and 
the last case (5). In fact in all the cases (l)-(3) the aspect ratio is very small (a„ <K 1), and so the 
in-plane interaction is stronger than the in-wall interaction, while the opposite situation arises in 
case (5), since a„ » 1 . For this reason we look at case (4) as the main critical case and we expect 
to get a better insight in the problem from its analysis. 

7.1. Heuristics for the scaling of the discrete energy 

Proceeding exactly as in the previous cases leads to the same requirement as in (36), where 
this time a n ~ 1 . Rewriting the condition (36) for an equidistant distribution of dislocations we 
have 

j, n n-k j, n 

i ~ y y v( a „k) — — y 

n z crha n j—[ j-J ncrha n 

Hence, since Yj™=\ < °°' tne boundedness of the energy reduces to the following condition 
on the coefficients: 

— \- ~ 1. (42) 
no~ha n 

From this bound we obtain the following expressions for the aspect ratio a n and the pile-up 
length {„: 

ncrh cr 

We note that the scaling regime can be rephrased as /3„ ~ 1 . Therefore we can rescale the energy 
(9) as in (38), i.e., 

4 4) (x) := § J] f V (np n ( Xj+k - Xj j) + - J] x jt (44) 

k=l 7=0 " 7=0 

where now /3„ ~ 1 (note that (43) suggests 0^ instead of /3„, but in this regime /3„ ~ fy. 

7.2. Continuum limit: derivation of the internal stress 

The derivation of the continuum energy in this regime is quite different from the one outlined 
in Sections 4-6. This can be explained by considering the discrete energy (9). While for a„ « 1 
(namely in the scaling regimes (l)-(3)) the sum over k gives rise to an integral term, in addition 
to the integral generated by the sum over j, for a n ~ 1 this is no longer the case, so a different 
approach has to be followed. We mention here the idea only briefly (setting /?„ = 1 for simplicity) 
and refer to [9] for the complete proof. 

The idea in this case is to view the positions of the walls x, as the deformed positions from 
an initial equispaced wall configuration of n + 1 walls in (0, 1) via a deformation map £„, i.e., 
Xi = ^n(^)- The map g„ can be extended in the whole range (0, 1) as the continuum affine 
interpolation of X\, . . . , x n (see Figure 1 1). 
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Figure 1 1 : Interpolation of the positions of the walls. 



In this way we can formally read the argument of V in (44) as a quotient ratio of f„, i.e., 
Hence for the discrete energy (44) we have 

*=1 { J=0 ) 7=0 

which for n — > oo converges to the continuum functional 

£ (4) ©:=X V(k?(s))ds+ {(s)ds, (46) 
frfJo Jo 

where £ is the limit of the interpolations and the integrals are the limits of the Riemann sums in 
(45). Since we are interested in a continuum model in terms of the dislocation density rather than 
the limit positions it remains to change variables in the energy (46), using the relation 

p(t(s)) := t!~ v (47) 

which is the continuum version of the discrete relation (5). The relation (47) entails ds = p(x) dx. 
Changing variables in (46) leads to a continuum energy depending on the dislocation density, 
namely 

£(4)(p)= X Ves (ptx)) p(x)dx+ L xp(x)dx ' (48) 
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where V eff (f) := ££=i V(JtO. for every t e M. 

The dislocation density p minimising the energy (48) is the solution of the Euler-Lagrange 
equation associated to the dimensionless functional £ (4) , i.e., 

-^ ff ^) + 1=0 ' ^ 
where <p e s(t) := k<p(kt), for every f e R. Note that the infinite sum that appears in the 
definition of the continuum equilibrium equation (49) results form having taken into account all 
the dislocation interactions in the discrete model. More precisely, the term k = 1 in the sum 
corresponds to the contribution of the interactions among nearest neighbours, the term k = 2 to 
next-to-nearest neighbour interactions and so on. Moreover, whereas in the previous cases the 
dislocation walls were sufficiently close to regard them as a continuum density, here the discrete 
interactions prevail. 

Equation (49) is again the mesoscopic equilibrium equation o^-cr = in its non-dimensional 
form. Hence, the dimensionless internal stress obtained for this rescaling is 

(4) d x p , l\ 
The dimensional version is (14). 

Remark 7.1. We show here how (49) can be heuristically derived starting from the original, 
dimensional equilibrium equations (3). We will focus on a general k-th term in the definition 

of <p eS . 

For a dislocation located at a point x, we denote by di the average distance to its k-th neigh- 
bours, i.e., dk ~ |. The actual distance df(x) between the dislocation at x and its k-th neighbour 
on the right and on the left are given by the following corrections of dk, namely 

k ( i dp\t\ k\ 

k I 1 dp\( 1 k 
d ^-p + \-p-^J\-2-p 




Figure 12: ^-interactions. 

Therefore, the distance between its k-th neighbours is 

Ad 



-^o- x p. 



According to the previous relations, the force exerted on the dislocation at x by the k-th neighbour 
dislocation on the right and on the left are, respectively, 

I - — u-l — 1 4- , 

2h r v /i 



K 



4)-fWr)-sWTrW 



h 
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The net force exerted on the dislocation at x by its £-th neighbours is, therefore, 



K K , K k n ,tk\ 

which is exactly the k-th term in (14). 

7.3. Comparison Discrete vs Continuum 

In this section we show the agreement between the solution of the continuum equation (49) 
and the minimiser of the rescaled discrete energy (44), for f3„ ~ 1 . The agreement is shown in 
Figure 13; also in this case we can notice the occurrence of boundary layers at the left end of the 
pile-up region. 
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Figure 13: Comparison of the discrete and continuous pile-ups for and i?' 4 ', with n = 150 and fl„ = 1. 



8. Supercritical regime a n » 1: derivation of the upscaled internal stress 

In this scaling regime the in-plane dislocation distance is much larger than the in-wall dis- 
tance, and therefore the in-plane interaction is quite weak. In this case we can prove that truncat- 
ing the in-plane interactions to the first neighbours (but keeping the walls) leads to the right limit 
model, unlike the previous cases, where all interactions had to be accounted for. 

8.1. Heuristics for the scaling of the discrete energy 

We first note that the energy (9) is always larger than its truncation to first neighbours, i.e., 

En(X) ~ ^te n % V ( na " {Xj+1 ~ Xj) ) + ~n Xj > (50) 
' 23 



regardless of the scaling regime. In the case a„ » 1 the bound (50) is optimal, meaning that 
the energy E„ and its first order truncation give rise to the same continuum model (see [9]). This 
is in contrast with what happens in all the previous cases (l)-(4), where the truncation of the 
interactions would produce a completely different (and wrong) result (see also the discussion 
in [22]). 

To obtain the energy scaling we use once more the equispaced wall configuration x, = ^ as a 
test configuration. Then the bound on the energy (actually the bound on the energy truncated to 
the first neighbours) reduces to 

A V(a„) ~ 1. (51) 



V(s) ~ -se- lns as s -> oo. (52) 



ncrha n 

Since a„ » 1, we can substitute for V its asymptotic behaviour at infinity; by (7) we have 

2 
n 

Therefore, using the previous relation, the bound (51) becomes 

— -e- 2 ™»~l. 
ncrh 7i 

In terms of a„, and using (20), the previous relation entails 

■?4M§«f)-sM?)- «> 

Therefore, the energy scaling for or„ » 1 can be rephrased in terms of f3 n as /?„ » 1, and the 
corresponding energy scaling in the regime f3 n » 1 is 

For what follows it is convenient to introduce a new parameter y„ : = ^ log^^- j (note that this 
scaling regime corresponds to y n » 1) and to rewrite (54) in terms of y„ as follows: 

S.2. Continuum limit: derivation of the internal stress 

As in Section 7 we introduce the interpolation f„ of the positions of the walls. Heuristically 
we expect the term k = 1 in the sum in (55) to be dominant, and hence we disregard the other 
terms in our formal derivation of the limit energy (for the rigorous proof we refer to [9]). We 
have, formally, 

n-l 



' J=0 J=0 
2iry„ nl r\ 

2 y„ Jo Jo 
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By (52) we have that, for s e (0, 1), 



2 y„ 



V(y n f'(s))~?(s)e 



2*(y„-fCs)) 



which is finite (actually zero) only if £'(s) > 1. In terms of the density p, using (47) the limit 
energy is the dimensionless continuum functional E (5> defined as: 



£ (5) (P) 



Jo 



xp(x)dx 



+oo 



ifp < 1, 
otherwise. 



(56) 



This limit functional is degenerate: it only is finite if the dislocations are sufficiently far apart. 
The optimal density in this case is the constant density 1 from x — to x = 1, and density zero 
for x > 1 (see Figure 14 below). 

8.3. Comparison Discrete vs Continuum 

In this section we show the agreement between the minimiser of the discrete energy (54), for 
/3„ » 1 and for large «, and the minimiser of the continuum energy (56). 

In Figure 14 we have plotted the piecewise constant optimal continuum density and the op- 
timal discrete density for n = 150 and j3„ = 10 8 : we see that the discrete density approaches the 
value 1. 
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Figure 14: Optimal densities for E„ ' and with n = 150 and /?, 
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9. Comparison with previous models 

In this section we show how the expressions for the internal stress obtained in this paper, 
starting from an idealised discrete model, relate to some well-known models in the engineer- 
ing literature. Our results offer therefore a unifying approach for understanding these existing 
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models: They can all be derived by upscaling the same discrete model, but under different as- 
sumptions on the local arrangements of the dislocations, i.e., for different values of the aspect 
ratio /?„. 

First of all we list the dimensional internal stresses obtained in the scaling regimes (l)-(4) 
(we omit (5) since the equilibrium equation is too degenerate to provide a useful expression for 
the internal stress), i.e., 



(1) Subcritical regime B n <sc - : cr Ax) = —r I log \x - y\dyp{y) dy; 

n lnt 7T 2 J,, 

(2) First critical regime /?„ ~ ^ : af^(x) = -K j V (^^) d f p(y) dy; 



1 (3) Kh „ 

(3) Intermediate regime - «: B n <sc 1 : cr: Ax) = -—dip; 

n mt 3n 

(4) Second critical regime p n ~ 1 : ^in^ 1 ) = ^eff d *P- 



The back-stress cr. J coincides with the one derived by Eshelby, Frank, and Nabarro [7] and 
Head and Louat [16]. The starting point of both works is a system of discrete equilibrium equa- 
tions for n dislocations in one slip plane, which is different from the walls of dislocations that we 
consider. This translates into a system of discrete equilibrium equations of the form (3), with the 
single-slip-plane interaction potential \jj(f) = l/(n 2 t) instead of the the potential tp. We note that 
the interaction potential 4> is similar to <f for small values of f, but differs for larger values. 

This is consistent with the fact that the internal stress ovj has been obtained in the scaling 
regime f3 n <sc 1/n, where all the dislocation walls are confined in a region that is small compared 
to the slip planes spacing h. Hence the in-plane interactions between dislocations are much 
stronger than the in-wall interactions. An approximation by a single-slip-plane setup is therefore 
appropriate. A comparison between the two discrete models in terms of the equilibrium densities 
is shown in Figure 7. 

The back-stress <r\^ coincides, at least qualitatively, with the one proposed by Evers, Brekel- 
mans and Geers in [8]. In the quoted paper the authors derive phenomenologically an expression 
for the internal stress that is, like ours, the gradient of the dislocation density, up to a multiplica- 
tive constant (representing, in their case, the range of interaction of a dislocation). 

In the special case /?„ = cn~ 1/2 for some constant c, the simple ordinary differential equation 
(41) characterising the equilibrium dislocation density in this regime has also been obtained 
(following a different method) in the recent paper [14], starting from the discrete system (3). 

The other internal stresses cr®, cr^ and cr^ were not obtained so far. 

The expression for the internal stress proposed by Groma, Csikor, and Zaiser [13] is a special 
case. It has been derived under the assumption that the distance to a nearest neighbouring dislo- 
cation is independent of the direction in which it is found. In our formulation this corresponds 
to the second critical case (4), namely /3„ ~ 1 . Therefore it is interesting to compare the internal 
stress of [13] to cr^ above. 
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As it turns out, the internal stress of [13], which up to constants reads cr ( mt ' = -d x p/p, can 
be formally obtained from o"S| by making two approximations. The first consists in truncating 
the number of interacting dislocations to the nearest neighbours. As described above, in the 
discrete model of this paper we take into account all interactions. This can also be recognized in 
o~- t , which is defined in terms of (fes, which in turn is a sum whose k-th term corresponds to the 
interactions between k-th neighbours. Hence, truncating the interactions to the first neighbours 
is equivalent to replacing the effective potential <p e s with <p. This first simplification reduces the 
internal stress erf 4 ,* above to 

int 

^ = T^*'Wl d * p - (57) 

The second approximation that leads to the internal stress proposed in [13] is to substitute 
the force <p(t) with its first-order Taylor-Laurent expansion close to zero, which coincides with 
the single-slip plane force if/{t) = \ used in [7, 16]. By using iff instead of <p in (57) we obtain 
exactly (up to a constant) crf^ cz) = -dip I p. 

The comparison between the interactions stresses (p e s, ip and if/ is illustrated in Figure 15. 




Figure 15: Interaction potentials tp e( f, ip and ijj in a log-log plot. 

From the figure it is clear that ip is a good approximation of <^ e ff only for large values of its 
argument. Therefore the first approximation (57) is justified when hp is small. This is exactly the 
case, though, when the second approximation is not allowed, since ijj and ip are close only near 
zero (corresponding to the opposite case of large hp). 

Therefore, although this derivation can formally be made, it can not be made rigorous, since 
the two approximations are mutually incompatible. More precisely, if first truncating to nearest 
neighbours (replacing <p e s by ip) and then Taylor-expanding ip were a consistent combination, 
then the Taylor expansion of ip^ and of ip would be similar. However, this is not the case, since 
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close to zero <p(t) ~ ip(t) = 1 /jr 



■ 2 t, while 2 



i r° 



as t — > 0. 



(58) 



Close to zero, therefore, many neighbours are relevant, and nearest-neighbour truncation results 
in a large error. 

The idea of deriving a continuum dislocation model from a discrete model of dislocation 
walls in a pile-up is also present in the work of Mesarovic and collaborators [2, 19]. Their 
approach consists in performing a two-step upscaling from the discrete dislocation model: first 
the dislocations are smeared out in the slip plane, but the discreteness is kept inside the wall; then 
the semi-discrete model obtained in the first step is upscaled in the vertical direction. We stress 
that there is no physical reason for performing the upscaling in the two directions separately. In 
fact, as the authors of the quoted papers correctly point out, the error between the predictions 
of the continuum model they obtain and the discrete model they start with (they refer to it as 
"coarsening error") is significant (and this was already observed in [22]). Therefore, they need 
to correct ad hoc the resulting continuum model in order to match the predictions of the discrete 
model. 

The continuum model we obtain, instead, agrees perfectly with the discrete model we consid- 
ered, and therefore we do not have to add any artificial terms to guarantee the matching between 
the two. 

10. Comments 

In this section we collect a number of interesting comments on the results presented in the 
paper. 

Transition between the different limit models. The transition between the regimes (l)-(5) is 
continuous in terms of the energy. For example, the integrand in the first term of the continuum 
energy £ (2) in (33) contains the term c V(cs), which for c large (corresponding to a transition 
from regime (2) to regime (3)) converges to (J V) 6, which is indeed the term appearing in the 
energy £ (3) in (40). If instead c is small (corresponding to the transition from regime (2) to 
regime (1)), then the logarithmic singularity of V appears, leading to the energy (24). The other 
transitions can be explained analogously. 

Mechanical interpretation of/3„. The dimensionless parameter j3„ measures the elastic prop- 
erties of the medium (described by K) in comparison with the strength of the pile-up driving force 
cr. Large fi n , therefore, corresponds to weak forcing, and small j3„ to strong forcing. Stronger 



2 This follows from the two inequalities (since s h» S(p(s) is decreasing in (0, oo) and s<p(s) = V) 




and 
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forcing pushes the dislocation walls closer to each other; when n/3 n — > 0, the forcing pushes the 
walls so close to each other that their distance is always smaller than h, so that V is only sampled 
close to the logarithmic singularity. 

Similarly, when /?„ — > oo, the forcing is so weak that the distance between the walls falls 
in the exponential tails of V. The intermediate regime is characterized by dislocation walls that 
span the range from 'smaller than K to 'larger than h\ 

Length scales. In each of the scaling regimes (l)-(5) it is possible to express the length of 
the pile-up region in terms of the parameters K, cr, h and the number n of dislocations. As one 
might intuitively expect, the length of the pile-up region will increase when going from (1) to 
(5). This is consistent with the previous observation that small /3„ corresponds to strong forcing 
and to confining the dislocations in a small domain, while for large f3 n the forcing is weak, and 
consequently dislocations can spread out over a larger region. 

More precisely, in the scaling regime (1) we found in (19) that the length of the pile-up scales 
like so in particular it is independent of the wall spacing h. This is not surprising, since in 
this case only the in-plane interactions play a role, so the walls could be equivalently replaced 

Knh 



by individual dislocations. The length of the pile-up in the cases (2), (3) and (4) scales as 
The case (5) exhibits an even different length scale, namely nh log ( 

In terms of /3„ the length of the pile-up region in each scaling regime is given by = n 2 0^h, 

4 2 ~ 4) = nj3 n h and £^ = g log (^), respectively. 

We can explain the transition between the length scales {fp - ffi in the following way. We 
notice that 

4»<e 4) Pn<-- 

n 

Therefore the transition between ffl and ^ 2_4) happens exactly in the first critical regime (case 



(2)), where f3 n ~ \. Concerning the transition between 4) and ( { n \ although ^ 5) < 6„ 4) al- 
ways, 4i 5) is acceptable only when /3 n is at least order one (since otherwisi 
And this condition corresponds to the second critical regime (case (4)). 



More general dislocations arrangements. A natural extension of this work is to apply our 
rigorous upscaling procedure to more general dislocation arrangements. A first direction is the 
study of random walls, namely perfectly straight (possibly finite) dislocation walls where the 
spacing between dislocations is not constant. As pointed out in [14] the problem in this more 
general case is genuinely two-dimensional, and hence much more difficult to treat. We refer also 
to [25] for a related discussion on random walls of dislocations. 



11. Conclusion 

This paper unravels the mechanical response of a system of walls of parallel edge dislocations 
that move along equidistant parallel slip planes. By implementing a rigorous mathematical limit 
procedure, we have identified five different parameter regimes, as the number n of walls tends to 
infinity. These regimes are characterized by the asymptotic behaviour of a single dimensionless 
parameter. 
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For each of these regimes we have identified the limiting internal stress that is generated by 
the dislocation density, and the structure of a pile-up of walls against a hard obstacle. For two of 
the regimes the expressions that we obtain were known in the literature, and our results provide 
new insight by clearly delineating the conditions under which these expressions are valid. The 
three other cases are new, and the corresponding behaviour has not been studied before. 

Although the analysed wall configuration is highly simplified, the rigorous nature of this 
work implies that these results can be considered as a benchmark: any (possibly more general) 
model that describes the behaviour of large numbers of dislocations should reproduce at least the 
behaviour given by the results of this paper when applied to the corresponding idealised situation. 
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